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Abstract 

We comment on the first indication of geo-neutrino events from KamLAND 
and on the prospects for understanding Earth energetics. Practically all mod- 
els of terrestrial heat production are consistent with data within the presently 
limited statistics, the fully radiogenic model being closer to the observed value 
(~ 9 geo-events). In a few years KamLAND should collect sufficient data for 
a clear evidence of geo- neutrinos, however discrimination among models re- 
quires a detector with the class and size of KamLAND far away from nuclear 
reactors. We also remark that the event ratio from Thorium and Uranium 
decay chains is well fixed N(Th)/N(U) ~ 0.25, a constraint that can be useful 
for determining neutrino oscillation parameters. We show that a full spectral 
analysis, including this constraint, further reduces the oscillation parameter 
space compared to an analysis with an energy threshold E v i s > 2.6 MeV. 



I. INTRODUCTION 

Recently KamLAND presented the first results [1] on the search for oscillation of v e 
emitted from distant power reactors. Electron antineutrinos are detected by means of inverse 
beta decay, 

v e +p -> e + +n - 1.80 MeV , (1) 

by looking at the prompt energy deposited by the positron (E vis = 2m e + E^m, where the 
kinetic energy of the positron is E kin = E p — 1.80 Me V ) accompanied by the signal of 
the neutron from n + p — ► d + 7. With an exposure of 162 ton-yr a clear deficit has been 
observed, however various combinations of oscillation parameters describe well the shape of 
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the positron spectrum. The best fit value reported in [1], including the geo-neutrino fluxes 
as free parameters, corresponds to sin 2 29 = 0.91 and Am 2 = 6.9 • 10" 5 eV 2 . 

Terrestrial antineutrinos, emitted by the /3-decay of the progenies of 238 U and 232 Th in the 
Earth's interior, contribute to the low energy part of the detected signal, the maximal E vis 
being 2.48 and 1.46 MeV respectively. From a fit to the experimental data the KamLAND 
collaboration reported 4 events associated to 238 U and 5 to 232 Th [1]. These numbers provide 
a direct insight on the radiogenic component of the terrestrial heat. In this paper we 
comment on the implications of this first result and on the prospects which it discloses 
for understanding the energetics of the Earth. We also discuss the constraints provided by 
geo-neutrinos for precise determinations of the neutrino oscillation parameters. 



Given the cross section of (1) and the antineutrino spectra one can immediately derive 
the relationship between the number of events and the antineutrino fluxes (see Appendix): 



where event numbers N are calculated for an exposure of 10 32 protons-yr, P is the averaged 
survival probability and e are the detection efficiencies (from the values quoted in [1], we 
get P = 1 - 1/2 sin 2 29 ^ 0.55 and e = 78.3% for both U and Th). $(X) are the produced 
fluxes in units of 10 6 cm~ 2 s _1 , i.e. the fluxes which one should observe in the absence of 
oscillation: 



where R is the detector location, p is the density, Cx, Tx and mx are the concentration, 
lifetime and atomic mass of element X and nx is the number of antineutrinos emitted per 
decay chain. The integration is performed over the Earth volume V®. 

The radiogenic contribution to the terrestrial heat is not quantitatively understood. In 
ref. [2] three representative models have been considered: 

a) a naive chondritic model, where one assumes for the Earth mass ratios typical of 
carbonaceous chondrites [3]: [Th\/[U\ = 3.8, [K\/[U\ = 7 • 10 4 and [U\/[Si\ = 7.3 • 10~ 8 . In 
this model the radiogenic heat production rate is about 30 TW, originating mainly from K 
decays. 

b) The Bulk Silicate Earth (BSE) model, which provides a description of geological 
evidence coherent with geochemical information and accounts for a radiogenic production of 
about 20 TW. In this model one has: [Th\/[U\ = 3.8, [K\/[U\ = 10 4 and [U]/[Si] = 9.4-10~ 8 . 

c) A fully radiogenic model, where the abundances of Th, U and K are rescaled with 
respect to b) so as to account for the full terrestrial heat flow of 40 TW. 

Uranium mass in the crust M crust (U) = 0.4 • 10 17 kg has been fixed. By taking 
M(Si)/Mq = 0.15 all other masses in the crust and in the mantle are obtained from the 
above ratios. Uniform distributions within the crust and the mantle are assumed. For each 



II. KAMLAND AND TERRESTRIAL HEAT SOURCES 



N(U) = 13.2 • P ■ e(U) ■ $(U) 
N(Th) = 4.0 • P ■ e{Th) ■ $(Th) 



(2) 
(3) 




(4) 
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model, from the fluxes of ref. [2] and eqs. (2,3) we obtain the number of events expected in 
the first exposure of KamLAND, see Table I. 

In view of the limited statistics, it is useful to consider the sum of terrestrial events 
N(U + Th). The measured value is essentially obtained from a total of C = 32 counts with 
E vis < 2.6 Me V, after subtraction of reactor events R and background B: N(U + Th) = 
C — R — B. The statistical fluctuation is thus of order \/C = 5.7. Within this uncertainty 
all models are consistent with data, the fully radiogenic model being closer to the value 
reported in [1]. 

We remark that the ratio N(Th)/N(U) is a significant indicator. In fact, the separate 
event numbers depend on the amount of radioactive materials, on their distribution inside 
the Earth, on the antineutrinos survival probability and on the detection efficiency. On the 
other hand, if one assumes an approximately uniform mass ratio [77i]/[C/] inside the Earth, 
the event ratio does not depend on material distribution and on the survival probability, as 
it is clear from eq. (14). Assuming e(Th) = e(U), one has 

N(Th)/N(U) = 0M5[Th]/[U] . (5) 

For most models of terrestrial composition one has [77i]/[C/] = 3.8, giving 

N(Th)/N(U) = 0.25 . (6) 

On the other hand, by considering both N(Th) and N(U) as independent parameters 
KamLAND obtains N(Th)/N(U) ~ 1. If confirmed with higher statistics, this would imply 
[Th}/[U] w 16, quite an unexpected value. However, a model with [Th}/[U] = 16, [if]/[C/] = 
10 4 and M crust+mant i e (U) = 0.8 • 10 17 kg would provide the full observed heat flow, the main 
source being Th at 28 TW. The model predicts about 5 events in KamLAND, half from Th 
and half from U decays. 

In order to discuss the achievable improvements, we collect in Table II the predictions 
of several models, normalized to 10 32 protons-yr v ) and assuming 100% detection efficiency. 
We present the results following from [2] , together with predictions obtained from the fluxes 
estimated in [6] and some estimates from [5], rescaled for a 0.55 survival probability (model 
lib of [5] assumes that heat production is fully sustained by U and Th, omitting any con- 
tribution from Potassium). From the various models, one estimates Th + U events in the 
range 24-83 respectively. 

At a site with negligible reactor flux N(Th + U) could thus be measured with an accuracy 
of about 20-10% and the different models could be clearly discriminated. On the other hand, 
by rescaling the present KamLAND data, one expects that counts with E vis < 2.6 MeV will 
be C ~ C ■ 7.19/0.78 = 295. This implies statistical fluctuations of about ±17 events, which 
possibly will allow for a clear evidence of geo-neutrinos, however they are too large for model 
discrimination. 

All this calls for a detector with the class and the size of KamLAND, far away from 
nuclear reactors. We note that BOREXINO [7] will provide additional and complementary 



KamLAND will presumably obtain this exposure within two years. We remind that KamLAND 
present fiducial mass is 408 tons, out of the total of 1000 tons of mineral oil. 
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information in the future. Its target mass is about 300 tons and the reactor background 
corresponds to about 7 events per year below 2.6 MeV, thus providing a better signal to 
background ratio [5,6]. 



III. GEO-NEUTRINOS AND OSCILLATION PARAMETERS 



When KamLAND data at E vis > 2.6MeV are combined with solar and Chooz data, 
the solution to the solar neutrino problems is basically split in two near regions, called 
LMA-I and LMA-II [8-14]. The first region contains the global best fit point, correspond- 
ing to (Am 2 /10~ 5 eV 2 , sin 2 6) = (7.3,0.315), whereas the second one is centered around 
(15.4,0.300) [8]. A relevant question is thus if geo-neutrinos can be of some help for dis- 
criminating between the two solutions. 

As previously remarked, although the total amounts of U and Th inside Earth are not 
well determined, the ratio of their abundances is rather constrained. Estimates for the solar 
system yield [77i]/[C/] = (3.7 — 3.9) [4], estimates for the primitive mantle are in the range 
(3.6 — 4), measurements of the upper continental crust give (3.8 — 4.2), estimates of the bulk 
continental crust are in the range (3.8 — 5) [15,16]. By assuming [77i]/[C/] = 3.8 ±0.7, from 
eq. (5) we get for the ratio of geo-events: 



We remark that this constraint, which has been derived by assuming an approximately 
uniform distribution of [T7i]/[[/] and equal (distance averaged) survival probabibilities, has 
actually a larger validity. 

Concerning the effect of regional [T7i]/[[/] variations, from [6] we derive that r is changed 
by less than 2% when the detector is placed at Kamioka, or Gran Sasso, or Tibet (on the 
top of a very thick continental crust) or at the Hawaii (sitting on the thin, U- and T/i-poor 
oceanic crust). Coming to the effect of local variations, by assuming that within 100 km 
from the detector the Uranium abundance is double, [Th]/[U]=2 , one gets r = 0.22, whereas 
if its is halved, [T7i]/[[/] = 8, one finds r = 0.28. Neutrino oscillations clearly do not affect 
eq. (7) if the oscillation lenghts for both U and Th neutrinos are both very short or very 
long in comparison with some typical Earth dimension. We have checked that the effect of 
finite oscillation lengths does not change r by more than 2In conclusion, all these effects are 
well within the estimated 20% uncertainty on r. 



In order to see the implications of this constraint, let us first divide the KamLAND signal 
below 2.6 MeV in two regions: a) 0.9 < E vis (MeV) < 1.75 corresponding to the first two 
bins of [1] and b) 1.75 < E V i S (MeV) < 2.60, corresponding to the next two bins. 

All T/i-events are contained in region a), whereas a fraction s of U events are in a) and 
in (1 — s) are in b). We find s = 0.6. The number of geo-events G a ^ in each region is thus: 



r = N(Th)/N(U) = 0.25 ± 0.05 



(7) 



A. A sum rule 



G a = (r + s)N(U) 



G b = (1 - s)N(U) 



(8) 
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By eliminating N(U) from the two equations one obtains a sum rule: 

(l-s)G a -(r + s)G b = . (9) 

For each solution, we can extract G a ,b from KamLAND counts (7 0j &, after subtracting the 
estimated background B a h and the predicted reactor events R a ,b- One can then build the 
quantity S = (l — s)G a — (r + s)G b and check if it is consistent with zero. We remark that C a ^ b 
are essentially independent observables and the statistical fluctuations are not correlated. 
This procedure is shown in Table III for the best fit points of LMA-I and LMA-II 2 . The 
resulting values of S, calculated for r = 0.25 and s = 0.6, 

S(LMA-I) = 1.7 ±3.7 ; S{LMA — II) — 3.6 ± 3.7 (10) 

show that both solutions are consistent with the sum rule. The constraint (9), which is prac- 
tically unaffected if r is varied within its assumed ±20% uncertainty, will become relevant 
when more data are available. 



B. Full spectral analysis 

In addition to the algebraic approach described above, we performed a fit to the entire 
positron spectrum (E v i s > 0.9 MeV) including the geo-neutrino contribution with [T/i]/[£7] 
fixed at 3.8. The fitting function includes then the reactor fluxes of the 16 main contributing 
power plants, as well as the geo-neutrino spectrum of 238 U and 232 Th [18] convoluted with 
the KamLAND energy resolution of 7.5%/^/E{MeV). Only the 5 first bins have a non-zero 
geo-neutrino contribution. According to [1] we included 2.9 background events; since the 
exact background distribution has not been published, we added 2 of them into bin 1 and 0.3 
into bin 2, 3 and 4 respectively. We renormalised the no-oscillation spectrum to 86.9 events 
for E vis > 2.6 MeV in order to match the KamLAND integrated exposure. This leads to 
about 122 expected reactor events for E vis > 0.9 MeV in the absence of oscillations. It is 
worth noting that in addition to the error of the overall normalisation (5.6 events), the lack 
of knowledge of the individual running time of the reactors adds another systematic error 
that we do not take into account. The x 2 function is taken as in [11] which accounts for bins 
with low statistics. Fitting of the KamLAND spectrum leads to two main minima which we 
label low-LMA and high-LMA 3 . Both solutions remain stable when increasing the threshold 
to E vis > 2.6 MeV; the best fit values do shift only slightly when including the additional 
geo-neutrino information. The detailed results of the fit are given in Table IV and displayed 
in Figs. 1 and 2. The best fit is obtained for the low-LMA solution (x 2 = 6.0/14) and the 



2 As a general consideration, we use here the values obtained from global analysis [8] and omit 
uncertainties related to theoretical predictions. 

3 For clarity, we refer to LMA-I and LMA-II when discussing results of combined analysis of 
KamLAND (E vis > 2.6 MeV) + Chooz + solar data, we refer to low-LMA and high-LMA when 
considering the energy spectrum of KamLAND. 
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corresponding geo-neutrino contribution N(U + Th) is found to be 9.9 ±6.2. The high-LMA 
solution is allowed at a one sigma level with A% 2 = 2.2 and N(U + Th) is 8.4 ± 5.9. This 
result can be compared with 9 events obtained in [1]. 

We applied the algebraic method described in the previous section using the low-LMA 
and high-LMA solutions, and find Si ow -lma = 0.6 ± 3.7 and Shi g h-LMA = 4.8 ± 5.0. This 
shows that the sum rule can be applied to check the consistency of the data and the geo- 
neutrino predictions. 

Finally we compare the constraints on oscillation parameters obtained from the analysis 
with E vis > 2.6 MeV and E vis > 0.9 MeV. In the latter case we include geo-neutrinos with 
the ratio [Th]/[U] fixed at 3.8 and background as described above. We calculate the values of 
X 2 in the 3-dimensional parameter space [Am 2 , sin 2 29, N{U + Th)]. To obtain the 95% C.L. 
for the subspace of interest, we then project the volumes which satisfy \ 2 ~ Xmin < 5. 99 (joint 
estimation of 2 parameters) onto the [Am 2 , sin 2 29} plane [20]. Fig. 3 shows the 95% contour 
lines of the full spectral analysis together with the contours obtained for E vis > 2.6 MeV. 
The latter is in good agreement with [1], however we remark that our 90% C.L. contour is 
closer to their 95% C.L.. 

Even with the present limited statistics, the full spectral analysis further reduces the 
allowed oscillation parameter space compared to the analysis with E vis > 2.6 MeV. We 
expect that with an increased statistics in the future, a full spectral analysis including geo- 
neutrinos will provide a severe consistency check of the data and moreover can help to break 
the degeneracy among the solutions, in particular, if a high-LMA solutions [19] were realized 
in nature. 
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APPENDIX: FROM FLUX TO SIGNAL 



The number of events N(X) from the decay chain of element X — U, Th is: 

N(X) =N p tJ dE p e(E p )a(E p )^ rr \Ep) (11) 

where N p is the number of free protons in the target, t is the exposure time, e is the detection 
efficiency and (px Tr \E p ) is the differential flux of antineutrinos arriving into the detector: 

Jv & 4ir\R-r\ 2 r x m x 

where p is the density, Cx, Tx and mx are the concentration, lifetime and atomic mass 
of element X and n x is the number of antineutrinos emitted per decay chain. f x (E p ) is 
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the energy distribution of emitted antineutrinos, normalized to unity, and p is the survival 
probability for v e produced at r to reach the detector at R. 

In view of the values of the oscillation length one can average the survival probability 
over a short distance and bring out of the integral the term: 

P = ( p ) = 1 - i sin 2 20 . (13) 

In this way we are left with: 

N(X)=N p tP fdE p e(E p )a(E p )f x (E p ) [ Sr ^ . (14) 

J Jv<$ 47r|ix — r\ 2 T X m x 

The second integral is the produced flux of antineutrinos $(X) of eq.(4). Also one can 
assume the detection efficiency as approximately constant over the small 1 MeV) energy 
integration region. This leads to: 

N(X) = N p tPe(X)$(X) J dE p a(Ep)fx(E p ) (15) 

This integral is easily computed from the cross section given in Ref. [17] and the spectrum 
from Ref. [18]. 
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TABLES 



TABLE I. KamLAND results and theoretical predictions. Events are estimated from 
[2], for 1.39 • 10 31 protons-yr, 78.3% efficiency and 0.55 survival probability. 





N{Th) 


N(U) 


N(Th + U) 


KamLAND 


5 


4 


9 


Chondritic 


0.53 


2.05 


2.58 


BSE 


0.62 


2.45 


3.07 


Fully radiogenic 


1.03 


4.03 


5.06 



TABLE II. Estimated geo-neutrinos events for 10 32 protons-yr, 100% efficiency and 0.55 



survival probability. 


Model 


[ref] 


N(Th) 


N(U) 


N(Th + U) 


Chondritic 


[2] 


4.8 


18.9 


23.7 


BSE 


[2] 


5.7 


22.5 


28.2 


Fully Radiogenic 


[2] 


9.5 


37.0 


46.5 


from KamLAND data 


[1] 


45.9 


36.7 


82.6 


from 


[6]* 


5.5 


21.1 


26.6 


la 


[5] 


8.6 


32.5 


41.1 


lb 


[5] 


5.7 


21.6 


27.3 


lib 


[5] 


14 


54 


68 


* Values obtained from the fluxes of [6] and eqs.(2,3). 
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TABLE III. Geo-events expected for LMA-I and LMA-II and the [Th]/[U] con- 
straint. The best fit parameters (Am 2 , sin 2 29) are from the combined analysis of KamLAND 
(E V i S > 2.6 MeV), solar and Chooz data from [8]. Survival probabilities P a ,b of reactor antineu- 
trinos from [8], counts C a ^, estimates of no-oscillation reactor events R a ^{n.o.) and background 
B a b from [1], Reactor events are estimated from R = P ■ R(n.o). The extracted geo-events are 
G = C — R — B. Errors on the measured counts C a ^ correspond to statistical fluctuations. 

Solution LMA-I LMA-II 



Am 2 (W- 5 eV 2 ) 


7.3 


15.4 


sin 2 29 


0.863 


0.840 


Pa 


0.65 


0.60 


Pb 


0.50 


0.58 


R a (n.o.) 


10 


10 


R b (n.o.) 


27 


27 


Ra 


6.5 


6 


Rb 


13.5 


15.5 


C a 


17 ±4.12 


17 ±4.12 


c b 


15 ±3.87 


15 ±3.87 


B a 


3 


3 


B b 








G a = C a — B a — R a 


7.5 ±4.12 


8.0 ±4.12 


Gb = Cb — B — Rb 


1.5 ±3.87 


-0.5 ±3.87 


N(U + Th) = G a + G b 


9 ±5.7 


7.5 ±5.7 


S = 0.4G a - 0.85G b 


1.7 ±3.7 


3.6 ±3.7 



TABLE IV. Energy spectrum analysis with and without the geo-neutrinos con- 
straint. Results from the full spectrum with [T7i]/[£7] = 3.8 constraint are compared with those 



from E vis > 2.6 MeV. 



Range of fit 




> 0.9 MeV 


Evis 


> 2.6 MeV 


Solution 


low-LMA 


high-LMA 


low-LMA 


high-LMA 


Am 2 (10" 5 eV 2 ) 


6.8 


14.8 


6.8 


14.8 


sin 2 29 


0.91 


0.84 


0.92 


0.78 


N(U + Th) 


9.9 ± 6.2 


8.4 ± 5.9 






x 2 


6.0 


8.2 


5.1 


6.9 


Data points 


17 


17 


13 


13 


d.o.f. 


14 


14 


11 


11 



10 



FIGURES 



— ■— KAMLAND data 

-~k- Best fit 

Reactors^l 
Geo^r s :U+Th 
Geo-% : Th 

- -f- - Geo^l : U 




6 7 8 
Prompt Energy (MeV) 



FIG. 1. Best fit to KamLAND data including geo-neutrinos with [Th]/[U] = 3.8 
(low-LMA). The various contributions to the sum spectrum are shown as derived by the fit, 
which gives N(U + Th) = 9.9 ± 6.2. The dotted vertical line corresponds to 2.6 MeV. 




Prompt Energy (MeV) 



FIG. 2. Second best fit to KamLAND data including geo-neutrinos with 

[Th]/[U] = 3.8 (high-LMA). N(U + Th) corresponds to 8.4 ± 5.9. 
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95.0% C.L., E Prompt >2.6 MeV 

I I i I i I i I i I i I i I i I i I i 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

sin 2 (2 0) 



FIG. 3. Comparison of the 95% C.L. allowed areas with and without geo-neutrinos 
constraint. The light-gray region is allowed analyzing the spectrum with an energy threshold 
of 2.6 MeV, while the dark-gray region is allowed including the full spectrum and constraining to 
[Th]/[U] = 3.8. 
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